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We consider classical normal modes and non-interacting bosonic excitations in disordered systems. 
We emphasise generic aspects of such problems and parallels with disordered, non-interacting sys- 
tems of fermions, and discuss in particular the relevance for bosonic excitations of symmetry classes 
known in the fermionic context. We also stress important differences between bosonic and fermionic 
problems. One of these follows from the fact that ground state stability of a system requires all 
bosonic excitation energy levels to be positive, while stability in systems of non-interacting fermions 
is ensured by the exclusion principle, whatever the single-particle energies. As a consequence, simple 
models of uncorrelated disorder are less useful for bosonic systems than for fermionic ones, and it is 
generally important to study the excitation spectrum in conjunction with the problem of construct- 
ing a disorder-dependent ground state: we show how a mapping to an operator with chiral symmetry 
provides a useful tool for doing this. A second difference involves the distinction for bosonic systems 
between excitations which are Goldstone modes and those which are not. In the case of Goldstone 
modes we review established results illustrating the fact that disorder decouples from excitations in 
the low frequency limit, above a critical dimension d c , which in different circumstances takes the 
values d c — 2 and d c = 0. For bosonic excitations which are not Goldstone modes, we argue that 
an excitation density varying with frequency as p(to) oc to 4 is a universal feature in systems with 
ground states that depend on the disorder realisation. We illustrate our conclusions with extensive 
analytical and some numerical calculations for a variety of models in one dimension. 



PACS numbers: 73.20.Fz, 63.50.+X, 75.30.Ds 



I. INTRODUCTION 

Excitations in condensed matter systems with 
quenched disorder have been a subject of intense study 
during the last several decades. Historically, it has been 
fermionic excitations in random systems that have re- 
ceived most attention. The reason for this lies in part 
with the rapid development of experiments and theory 
involving mesoscopic conductors, where the effects of dis- 
order in phase-coherent electron systems have been stud- 
ied in great detail. 1 

It is, however, also of considerable interest to study 
random systems with bosonic excitations, and there is an 
extensive literature treating problems of this type, too. 
For instance, the propagation of phonons in glasses and 
of electromagnetic waves in media with random refractive 
index has long been a subject of active research;^ and 
trapping of light via scattering from disorder is a princi- 
ple on which random lasers are based^ Other examples 
of bosonic excitations in random systems include vibra- 
tions of pinned elastic structures such as charge density 
wavesj^ magnons in diluted antiferromagnets and spin 



velop a framework for treating bosonic excitations which 
we use to calculate the density of states and localization 
properties of one dimensional, randomly pinned systems. 
Some of our results have been presented in short form 
elsewhere^ 



glasses r^S^ 



,10.11,12 



and quasiparticles in supcrfluid liq- 



uid helium permeating a porous medium^ To some ex- 
tent, work on these problems has focussed on specific 
features of individual examples, and given less emphasis 
to a generic aspects than has been the case for disordered 
fermionic problems. 

In this paper we emphasise just these generic aspects. 
In particular, we examine the relationship between uni- 
versality classes identified for fermionic problems and 
models for bosonic excitations, as well as features that 
are specific to bosonic problems. In addition, we de- 



Our discussion is organised as follows. In Seed] we 
review the symmetry classes established within random 
matrix theory for disordered fermionic systems. We re- 
call in Scc. lIIII the general form for a quadratic bosonic 
Hamiltonian and the Bogoliubov transformation required 
to diagonalise it, and show how it is useful to introduce 
an auxilliary problem with structure similar to that in the 
chiral symmetry class. In Sec . 11 VI we discuss how various 
examples of systems with bosonic excitations fit into this 
general framework. Here we emphasise the distinction 
between excitations that are Goldstone modes and those 
that are not, and outline an established argument that 
leads in the latter case to the result p(u>) cx u> . Then, as 
an illustrative example, we consider the random field XY 
spin chain, using our formalism as the basis for a numeri- 
cal study of the discrete system, presented in Sec.|3 and 
giving an analytical treatment of the continuum limit in 
Sec . 1 VII recovering in both cases the behavior p(oj) oc u 4 . 
In Sec. lVIII we consider the random field Heisenberg spin 
chain and some related but simpler models, which are of 
interest because disorder enters them in a more general 
way than for the XY model. Finally, we summarise the 
relevant experimental situation in Sec lVIlil and end with 
concluding remarks in Sec. lIXI 
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II. SYMMETRIES AND DISORDER 

To provide a context for our discussion of bosonic 
systems, we begin by setting out the symmetry classes 
that are recognized within random matrix theory for 
fermionic systems Models of non-interacting quasi- 
particles play an important role in the study of fermionic 
excitations in disordered conductors, insulators and su- 
perconductors. It has long been appreciated that the 
properties of these models are controlled by the discrete 
symmetries of the single particle Hamiltonians. Orig- 
inally three symmetry classes for random Hamiltoni- 
ans were identified>i£*i& One class consists of random 
Hamiltonians that have time-reversal invariance but no 
Kramers degeneracy. In an appropriate basis the Hamil- 
tonian H is real, so that 

H = H*. (1) 

Random Hamiltonians which obey Eq. appear, for 
example, when studying disordered conductors without 
applied magnetic field. A second symmetry class consists 
of time-reversal invariant random Hamiltonians for par- 
ticles with half-integer spin and hence Kramers degener- 
acy. In this case the time-reveral operation includes spin 
inversion, and invariance requires 

H = (j 2 H*a 2 . (2) 

(Here and in the following, <jj for i = 1, 2 or 3 represent 
the conventional Pauli matrices, acting on a subspace 
identified by the context.) In order for Eq. @ to be 
different in an essential way from Eq. spin rotation 
invariance must be broken. Thus this case is of relevance 
for disordered conductors with spin-orbit coupling. A 
third symmetry class arises when the Hamiltonian has 
no discrete symmetries. Examples of these symmetry 
classes are provided by the three Wigner-Dyson random 
matrix ensembles iLHL 

More recently, it has been recognizedi2i22i2ii22 that 
there are seven additional classes of disordered fermionic 
Hamiltonians. These arise where there exists a special 
reference energy (taken to be zero in the following) for the 
system, and a symmetry operation which relates eigen- 
states in pairs. Three of these seven are referred to as 
chiral symmatery classes ^StS&SL Hamiltonians for these 
classes can be put into the form 

where Q is itself a matrix or an operator. They obey the 
symmetry condition 

H = -a 3 Ha 3 . (4) 

This ensures that energy levels of such Hamiltonians ap- 
pear in pairs ±E, since if ip is an eigenfunction with 
energy E, then a 3 ip is an eigenfunction with energy — E. 
The symmetry condition of Eq. , when combined with 



either Eq. QJ, or Eq. J5J, or neither, leads altogether 
to three chiral symmetry classes. Chiral Hamiltonians 
appear as tight-binding models with only off-diagonal 
disordes2ii2ii22i2£i2& and in the problem of classical dif- 
fusion in a random medium^SS 

The remaining four symmetry classes arise in the study 
of disordered superconductors with pairing treated in the 
mean field approximation. The Hamiltonian for such a 
problem has the Bogoliubov-de Gennes structure 

where the kinetic term h is Hcrmitian, while the gap func- 
tion A is antisymmetric. This structure leads to another 
defining symmetry condition, 

H = -o-iifVi . (6) 

As for chiral Hamiltonians, the symmetry displayed 
in Eq. © ensures that energy levels of Bogoliubov-de 
Gennes Hamiltonians appear in pairs, ±E. There are 
four symmetry classes of such Hamiltonians, according 
to whether or not the system has time-reversal and spin- 
rotation symmetry, making the total count ten. 

A consequence of the conditions specified in Eqns. (@J 
and © is that statistical properties of energy levels and, 
in spatially extended systems, the associated eigenfunc- 
tions, are quite different in these additional symmetry 
classes near zero energy, compared to properties far from 
zero energy, or in the Wigner-Dyson symmetry classes^ 

It is natural to ask whether this classification can be 
extended to problems involving non-interacting bosonic 
excitations or, equivalently, classical normal modes. At 
first sight, it might seem that quasiparticle statistics are 
unimportant in a non-interacting system. In one crucial 
respect, however, this is untrue, since stability of a sys- 
tem requires bosonic excitation energies to be positive, 
while for fermionic excitations stability is guaranteed by 
the Pauli exclusion principle, whatever the single-particle 
energy levels. This has two implications. First, energy 
zero emerges as a special point in the spectrum of bosonic 
systems, as it does for the additional fermionic symme- 
try classes discussed above. And second, one anticipates 
that matrix elements of bosonic Hamiltonians for ran- 
dom systems will have specific correlations, which ensure 
positivity of the spectrum. Thus, while the most general 
form for a quadratic bosonic Hamiltonian (see Eq. (JJJ be- 
low) is superficially similar to the Bogoliubov-de Gennes 
Hamiltonian, and while normal mode frequencies, like 
the eigenvalues of Eq. JSJ, appear in pairs ±w, matrix el- 
ements must satisfy constraints in order that frequencies 
are real. Such a requirement is in stark contrast with the 
assumptions of statistical independence used in the con- 
struction of random matrix ensembles for fermionic sys- 
tems. To stress the significance of this point, imagine a 
treatment of a disordered, interacting system which pro- 
ceeds in two stages, by first finding the ground state and 
then calculating excitation energies within a harmonic 
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approximation. The spirit of random matrix theory for 
fermionic systems is to divorce these two stages and ap- 
proach the second one phenomenologically, choosing sta- 
tistically independent matrix elements. By contrast, for 
bosonic systems it is clear that the two stages cannot 
be completely separated. Indeed, whereas for random 
matrix theory universal spectral properties follow largely 
from symmetry and are independent of the details of 
the matrix element distribution, we argue here that for 
bosonic excitations which are not Goldsone modes it is 
the requirements of stability and the ensuing correlations 
in the Hamiltonian for excitations that give rise to uni- 
versal spectral properties. 



III. BOSONIC HAMILTONIANS 

In this section we discuss the most general form for 
bosonic Hamiltonians and summarise the diagonalisation 
procedure, following standard lines (see, for example, 
Ref.Eflh. We also emphasise the distinction between os- 
cillator frequencies and stiffnesses. Finally, we show for 
oscillations about a stable ground state that it is natural 
to rewrite the Hamiltonian in terms of a chiral matrix. 



A. Stiffnesses and frequencies 

The most general bosonic Hamiltonian can be written 
in the equivalent forms 

1 N 

ff = r X! [MijPipj + K lj q i q j + 2C i3 qiPj] 
= j(p q)( c t K ) ( q 



Here, qi and pi are the coordinates and momenta of the 
oscillators, and a|,aj = (qi ±ipi) j\[2 are bosonic cre- 
ation and annihilation operators. The matrices M and 
K are real and symmetric, while C is an arbitrary real 
matrix. Equivalcntly, T is hermitian, while A is sym- 
metric. Physically, M is the inverse mass matrix of the 
oscillators, K is the matrix of spring constants, and cou- 
plings of the type represented by C occur, for example, 
in spin systems. It is convenient to define the 2N x 2N 
symmetric matrix 



H = 



M C 
C T K 



(8) 



We are interested in frequencies of oscillators decribed 
by Eq. (JJJ . From Hamilton's equations of motion it fol- 
lows that these frequencies are eigenvalues of the non- 
Hermitian matrix 



7 (i i iC T — iK . _ .. 



(9) 



They are real if the system is stable, in which case H 
is bounded from below, the eigenvalues of TC are posi- 
tive, and we can write Ti. = Q T Q with Q real. In thse 
terms, to find frequencies we need to diagonalize the ma- 
trix a 2 Q T Q, but its eigenvalues coincide with those of 
another matrix 



n = Qa 2 Q T , 



(10) 



which is explicitly Hermitian and, moreover, antisym- 
metric, so that frequencies come in pairs ±tt>i, with 
i = l...N. 

If f2 is interpreted as a random Hamiltonian, then ac- 
cording to the classification scheme discussed in Sec.lTTlit 
belongs to one of the Bogoliubov-de Gennes classes (more 
precisely, to class D, see Ref.ll^). While this is indeed an 
indication that random oscillators behave in many ways 
like random fermionic Hamiltonians from one of the addi- 
tional symmetry classes (having frequencies in pairs, with 
lu = as a special point in the spectrum), the identifica- 
tion of fl as class D operator is not by itself necessarily 
helpful since f£ does not have statistically independent 
matrix elements, but rather is constrained to have the 
form given in Eq. 1)1(1 [) . Instead, we shall see that a link 
to matrices wih chiral symmetry proves more useful. 

The computation of oscillator frequencies can equiva- 
lently be described as a Bogoliubov transformation for 
coordinates and momenta, specified by real matrices g 
which obey 



<72 = gcr 2 g 



(11) 



and tranform Ti to gTig T ■ Diagonalizing Ti using this 
transformation, the Hamiltonian H of Eq. (7J takes the 
form 



(12) 



In addition to the frequencies, obtained as eigenvalues 
of or by Bogoliubov transformation of H, the eigenval- 
ues of Ti. also have physical significance, for example as 
inverse static susceptibilities. We refer to them as stiff- 
nesses, denoting them by Ki, i = 1, . . .2N. In general, 
there is no simple relationship between stiffnesses and 
frequencies, but several special cases provide important 
exceptions, as follows. Consider first H. with C = and 
M = 1, representing oscillations of particles, all with 
unit mass, connected by springs with spring constants 
Kij. In this example, half of all the stiffnesses are equal 
to 1, while the other half are the eigenvalues Ki of the 
matrix K; the frequencies and stiffnesses are related by 



(13) 



A second special case arises for magnon excitations in 
weakly disordered ferromagnets, which have a Hamilto- 
nian of the form of Eq. (7J with M = K and C = 0. In 
that case, the stiffnesses are the eigenvalues of M and K, 
and come in identical pairs. The frequencies are simply 



UJ, 



±Ki 



(14) 
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A final and important special case occurs when one stiff- 
ness, say K±, is much smaller than all others. In this 
regime, an approximate relation exists between the small- 
est frequency and the smallest stiffness: 



(15) 



To derive this relation, a more detailed analysis of the 
structure of TL and £1 is required. Since TL is a real sym- 
metric matrix with positive eigenvalues, it can in general 
be represented as 



n = ua 2 u 7 



(16) 



where U is an orthogonal matrix and Ay = \5ij is the 
diagonal matrix whose eigenvalues are square roots of 
the eigenvalues of TL. The corresponding matrix f2 can 
be written as 



fi = MJ T a 2 UK. 



(17) 



Introducing an antisymmetric matrix A = U T a2U, we 
have 



— \i Aij Xj 



(18) 



Suppose initially that one of the stiffnesses vanishes, so 
that Ai =0. Then at least one frequency must vanish 
as well, since det SI = detTL. In fact, because $1 is an- 
tisymmetric and its frequencies come in opposite pairs, 
two frequencies are zero. The mathematical mechanism 
for this is clear from Eq. (|18fl . First, since flu = flu = 
for all i, one of the eigenvalues of Q is ui\ = 0, with 
an eigenvector ip^ = du. Second, all other eigenval- 
ues can be found by diagonalizing a smaller matrix, fi^- 
where 2 < i,j < 2N. But fi' is an odd-dimensional an- 
tisymmetric matrix, and therefore has at least one zero 
eigenvalue, u)2, with an associated normalised eigenvec- 

(2) (2) 

tor which we write as ip\ , where ip\' = 0. Now treat 
small non-zero Ai using perturbation theory about this 
limit, with Ai = e. The change in f2 is 



e (SuA^Xj + XiAjSij) (1 - S u S l3 ) . (19) 



Under this perturbation, the doubly degenerate eigen- 
value ui = of $1 splits into ±u>\, determined by diago- 
nalizing the 2x2 matrix 





Ei= 2 Mi\il>. 



(2) 



(2) 



Hence 



u>i = ±Ai 



2iV 



y] AuXjip. 



(2) 



(20) 



(21) 



barring an accidental vanishing of the matrix clement 

The usefulness of this result lies in the following. Con- 
sider a random system with bosonic excitations that 



are localized with a finite localization length at low fre- 
quency, and let the density of stiffnesses be d(n). We 
expect that each localization volume can be treated as 
an independent system, and that each low frequency ex- 
citation will have a frequency much smaller than that 
of other excitations in its own localization volume, and 
will therefore be associated with a single small stiffness. 
Applying Eq. (|15fl , the density of excitation frequencies 
p(uj) for small u> is 



p(oj) — d(LU 2 )(jJ . 



(22) 



One interesting check of these conclusions is provided 
by random matrix theory. Consider an ensemble of real 
random matrices Q, with size N 3> 1 and probability 
distribution P oc exp (— Q T Q). Let A be the typical 
magnitude of the eigenvalue o{H = Q T Q closest to zero. 
For 7i generated in this way, one finds d(n) cx k^ 1 / 2 at 
both k <C A and k 3> A. Computing the eigenvalue 
density of fl = Qa^Q 7 using superintegrals, one finds 
thatSl p{oj) is independent of u> for uj <C A, while p{ui) oc 
U)- 1 / 3 for lu A. It is clear that in the regime w -C A, 
Eq. 122|) is indeed applicable, while in the opposite regime 
it breaks down. 



B. Chiral Symmetry 

To study stiffnesses of bosonic oscillators with H = 
Q T Q, it is advantageous to introduce 7Y, an auxiliary 
matrix in which Q and Q T enter linearly. 



n = 



Q 
Q T 



(23) 



Obviously, the eigenvalues Ai of TL are square roots of 
the stiffnesses Ki of TL, and the matrix TL plays the role 
of the square root of the original bosonic Hamiltonian 
TL. On the other hand, the off-diagonal structure of TL 
is the defining feature of the chiral symmetry class, dis- 
cussed in Sec. |n] While the direct implications of this 
connection are limited, because the elements of Q are not 
independent random variables as would be the case in a 
random matrix ensemble, techniques originally developed 
for systems in this symmetry class will prove useful in our 
treatment of one-dimensional systems, as we describe in 

Sec. El and Sec. Ell 

To study the frequencies of the oscillators, as opposed 
to their stiffnesses, a second auxiliary matrix 



H' = 







Q 

<r 2 Q T 



(24) 



is helpful. We shall call matrices with this structure chi- 
ral bosonic matrices. The eigenvalue equation for TL' can 
be written in the form 







1> = 





0"2 



(25) 
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which of course inherits its structure from Hamilton's 
equations. 

The main difficulty in making use of these ideas is that, 
in practice, only ri is known initially and a method must 
be developed to find Q. Moreover, Q is defined only 
up to a left multiplication by an arbitrary orthogonal 
matrix. For the introduction of Q to be helpful, it will 
be important that it can be chosen to have a simple form, 
with, for example, only short range couplings. We shall 
show that this is indeed possible for a variety of problems. 



IV. BOSONIC EXCITATIONS: GENERAL 
ASPECTS 

It is a feature of models for disordered fermionic sys- 
tems from the chiral and Bogoliubov-de Gennes symme- 
try classes that their characteristic behavior appears only 
close to the reference energy, identified by the discrete 
symmetry of the Hamiltonian, while spectral properties 
at energies far from this are indistinguishable from those 
of the Wigner-Dyson classes. In a similar way, for bosonic 
excitations in disordered systems we expect it to be prop- 
erties at low frequency that that are of particular interest. 
Examples of such excitations can be divided into two cat- 
egories, according to whether or not they are Goldstone 
modes, associated with a broken continuous symmetry. 
In this section we summarize for each of these categories 
some of the previously established results. We also illus- 
trate the introduction of a chiral Hamiltonian ri for a 
one-dimensional phonon model. 



A. Goldstone Modes 

Important instances of Goldstone modes in disordered 
systems include acoustic phonons in glasses and alloys* 2 ^ 
and spin waves in isotropic, dilute ferromagnets and an- 
tiferromagnets, and isotropic spin glasses&LSiSiiSiiiiiS In 
disorder-free versions of these systems, low frequency 
excitations have long wavelength and are described by 
equations of motion that involve macroscopic properties 
of the system: density, elastic constants, magnetic sus- 
ceptibility and spin stiffness. Disorder introduces local 
fluctuations in the values of these quantities, but one 
expects excitations to couple only to fluctuations aver- 
aged over a volume with linear dimensions set by a wave- 
length. Because of this, randomness only weakly affects 
low-frequency Goldstone modes, especially in higher di- 
mensions for which the averaging is most effective. Such 
averaging is demonstrated by the low-frequency behavior 
of the excitation density p(u>), which above a critical di- 
mension d c varies with the same power of oj as in the 
disorder-free system. It is also shown by localization 
properties of excitations: in one and two-dimensional 
systems, the localization length £(u>) diverges as u> ap- 
proaches zero, while in higher-dimensional systems, all 
low-frequency states are extended. In this subsection we 



review the behavior of p(u>) and £(w) for phonons in al- 
loys and for spin waves in diluted antifcrromagnets and 
spin glasses, considering the effect of weak disorder in- 
cluded in the relevant equation of motion. We also use 
the Hamiltonian for phonons in a one-dimensional disor- 
dered system as an illustration of the general mapping to 
chiral models, introduced in Sec. IIII Bl 

Consider first a scalar version of a model for acoustic 
phonons, in which a mode with frequency u> has coordi- 
nate q(r) satisfying 

uj 2 q{v) = -c 2 V 2 <?(r) . (26) 

Suppose that the speed of sound, c has random fluctu- 
ations in space, about an average value Co, with only 
short-range correlations. Work on this and related prob- 
lems is reviewed in Ref. 0, and an early treatment of 
localization in this context was given in Ref. |2j. The es- 
sentials for our purposes are as follows. First, the fluc- 
tuations in c, averaged over a d-dimensional volume of 
size set by the wavelength in the disorder-free system, 
decrease compared to Co with decreasing frequency as 
uj d / 2 . Thus in this case the critical dimension is d c = 0, 
and for any d > d c the excitation density approaches the 
form found without disorder at low frequency. Moreover, 
from a calculation using the Born approximation, the 
Rayleigh scattering rate r _1 varies as t' 1 cx uu d+1 . (This 
dependence combines a factor of to 2 from the frequency- 
dependence of coupling to disorder, and a factor of to d ^ 1 
from the density of final states for scattering processes). 
In one dimension, disorder localizes with a localization 
length proportional to the mean free path, which here is 
cot, and so£ 

£(lu)(xlu- 2 . (27) 

In two dimensions it is a familiar feature of electronic 
systems that the localization length is exponentially large 
in kpl (where kp is the Fermi wavevector and / is the 
mean free path) : the equivalent parameter for the phonon 
problem is lot, so that 2 

£H cx exp([w / W ] 2 ) , (28) 

where loq is a disorder-strength dependent constant. 

These results contrast interestingly with those for an 
antiferromagnet which has randomness generated either 
by site dilution (taken small enough that the system is 
above the percolation threshold) or by substitution of 
impurity spins with a magnetic moment different from 
that of the host spinsAii In a discussion of the random 
antiferromagnet, it is useful to begin from the dispersion 
relation for spin waves in a two-sublattice ferrimagnet 
without disorder. At small wavevector k, this has the 
form 

lo 2 + (M a - M b )u) = c 2 k 2 , (29) 

where M a — Mb is proportional to the difference between 
the two sublattice magnetic moments, which are taken 
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to be oppositely aligned. Setting M a = Mb, one recov- 
ers the usual dispersion relation in an antifcrromagnet, 
with spin-wave speed c. For the random antiferromag- 
net, independent disorder on the two sublattices gener- 
ates random fluctuations in the local value of M a — Mb 
about a mean value of zero. Averaging these fluctuations 
over a volume of size set by the wavelength in the undi- 
luted system gives a random variable with an amplitude 
that scales as io d l 2 . Because of this, disorder appears in 
Eq. J2H1) via the term (M a -M b )uj, which scales as uj d / 2+1 . 
Above the critical dimension d c — 2, disorder is irrelevant 
in the sense that this term may be be neglected for small 
u) compared to to 2 . The value of the critical dimension 
is also apparent from a Born approximation calculation 
of the rate for scattering of spin waves by disorder, 11 
which yields t _1 cx ui d ~ x . Because of this, spin waves 
have a well-defined wavevector in the low frequency limit 
for d > d c , since lot — > oo as to — > 0. Applying this 
approach below the critical dimension, for d — 1 we see 
that Eq. determines a relation between the length- 
scale of an excitation, which we denote by 1/k, and its 
frequency: 



k d / 2 uo cx k 2 



at small to. This implies 

p{bj) CX LU~ 1/3 

and 

CMcx^- 2 / 3 . 



(30) 



(31) 



(32) 



We expect all these results for low-frequency behav- 
ior to be characteristic not only of disordered antifer- 
romagnets, but also of spin glasses, since the two sys- 
tems have in common the crucial feature of a magneti- 
zation density that is locally non-zero and random, but 
has zero average. A detailed treatment of excitations 
in spin glasses, however, is much more difficult than 
in weakly disordered antiferromagnets, because for spin 
glasses the ground state is generally disorder-dependent 
and unknown. Instead, the established approach is a 
hydrodynamic onef&ifi which leads to linearly dispersing 
modes with a speed determined by the macroscopic spin- 
stiffness and susceptibility. In the light of our discus- 
sion, we expect these hydrodynamic results to be correct 
for d > d c . By contrast, for one-dimensional systems, 
microscopic calculations are possible since frustration is 
absent and ground states can be determined. Results 
from computations^ in one-dimensional models of p(w) 
and £(w) coincide with Eqns. (|31|) and (|32|l . above. Sim- 
ilar calculations are also possible in higher dimensions 
for the Mattis model, which shares with one-dimensional 
models the features that frustration is absent and the 
ground state is known for each disorder configuration. In 
d = 2 these yield p(uj) cx lo\ log(cj)|, where the logarithm 
is characteristic of behavior at a critical dimension, and 
in d = 3 they give p(io) cx to 2 , in agreement with hy- 
drodynamic theory^ A discussion of excitations in the 



Mattis model, building on the methods described in this 
paper, will be presented elsewhere^ 

Spin waves have also been investigated 34 in the infinite- 
range Heisenberg spin glass: one of the main findings is a 
density of states that varies with frequency as p(oj) cx lo^ . 
We have not been able to make contact between this 
result and the approaches described here. 

As a next step, it is interesting to return to acous- 
tic phonons in disordered systems and use the one- 
dimensional version of this problem to illustrate some 
of the methods set out in Sec. IIHI In fact, in this con- 
text the mapping to a chiral problem was exploited in 
celebrated early work by DysonjS and also in calcula- 
tions by Zimani 2 ^ Consider a one dimensional chain of 
particles with masses m„ connected by nearest-neighbor 
springs with spring constants ki, where rrii and ki are 
random and positive. Let Pi and qi be the momentum 
and displacement of the z-th particle The Hamiltonian is 



Pi 

2m, 



-fa 



2 



(33) 



It is convenient to make the canonical transformation 
Pi -> ^/rrTiPi and ^/m~, giving 



El 

2 



(34) 



which is a particular case of the general bosonic Hamilto- 
nian Eq. (JZJ), in which M — 1 and C — 0, so that eigenfrc- 
quencies are related to stiffnesses n that are eigenvalues 
of the matrix K, as in Eq. (|13fl . Further discussion is 
easier in the continuum limit: replacing the index i with 
a continuous coordinate x, the Hamiltonian becomes 



H 



dx 



p 2 (x) k(x) 



+ 



( d 


q(x) 


)'] 


\ dx 


_y/m(x)_ 





(35) 



To find the stiffnesses, one must solve the eigenvalue 
equation Kq(x) = nq(x) with 



Kq{x) 



1 



d -k{x) d 



q{x) 



y/m(x) dx dx yjm(x) 



(36) 



Note that the operator K must be positive definite, 
in order that the c hain is stable. A nd in deed, defin- 
ing a(x) — y/k(x) and b(x) = l/y/m(x), it may be 
expressed as a square, in the form K = Q T Q with 
Q = a{x){d/ dx)b{x). Introducing a chiral Hamiltonian, 
as in Eq. (|23|l . we then have 





-b(x)£a(x) 



a{x)±b{x) 




1p = UJljj . 



(37) 



At the equivalent point in his treatment, Dyson dis- 
tinguishes between different possible choices for the form 
of disorder. In one case, termed type /, Eq. (|37J) is ef- 
fectively replaced by Eq. (|75|l below (with {V(x)) = 0), 
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leading to the singularity of Eq. (|61|l in p(ui). An al- 
ternative, termed type II, retains instead Eq. 137[) with 
only short-range correlations in a(x) and b(x), represent- 
ing disorder that couples only weakly at small k, because 
it is multiplied by spatial derivatives. This yields^ a 
constant p{u>) at small to, in agreement with the general 
arguments set out following Eq. J2HJ). ZimanS has given 
a detailed discussion in this context of the consequences 
of different types of disorder^ 



B. Non-Goldstone Low Energy Excitations 

Without Goldstone modes, the very existence of the 
low-lying excitations on which we have focussed our at- 
tention is not guaranteed. In fact, as seems first to have 
been appreciated in the context of atomic vibrations in 
glasses, 37 disorder itself may provide a route to a gapless 
spectrum. The essential ingredients are that the ground 
state should depend on the disorder realization, and that 
excitations at low frequency should be localized by disor- 
der. Then it is reasonable to consider excitations within 
each localization volume separately, and to expect disor- 
der configurations that support low frequency excitations 
to occur with finite probability. Roughly speaking, these 
excitations occur in regions where the ground state con- 
figuration is unusally sensitive to small changes in the 
disorder. In this section we summarise an established 
approach to this phenomenon of disorder-generated low 
frequency excitations, which concentrates on a single co- 
ordinate and its conjugate momentum. In subsequent 
sections we apply the formalism developed in this paper 
to study the phenomenon more generally. 

Following Il'in and Karpov^ consider a one- 
dimensional anharmonic oscillator with momentum p, 
mass m, coordinate q and potential U(q). The Hamil- 
tonian is 



H=2- + U(q) 
Zm 



(38) 



Choosing U (q) to be a smooth random function, we wish 
to study the frequency of small-amplitude oscillations 
about the absolute minimum in U (q) and, specifically, 
the probability distribution of this frequency. To give a 
more precise meaning to the notion of a smooth random 
function, we expand it in Taylor series 



U(q) = J2' 



q" 



(39) 



n=l 



where, to fix the zero of energy, we set U(0) = 0. We take 
the a n to be random with joint probability distribution 
P{a\,a 2 , . . .). We shall assume that P(a\, a 2 , . . .) is free 
of zeros and divergences, but its detailed form will not 
be important. 

Suppose that U(q) has a minimum at q = qo. In order 
to discuss excitations of the oscillator, we first Taylor 



expand U(q) around q — qo, writing 

■{q-q ) n (-go)' 



U{q)=Y,bn 



71=2 



(40) 



The probability distribution of the coefficients b n is re- 
lated to that for the a n by 



P(qo,b 2 ,b 3 ,. 



det 



d{ai,a 2 ,a 3l . . .) 



P(ai,a 2 , a 3 , . . .) . 



d(qo,b 2 ,b 3 , . . .) 

Evaluating the Jacobian, we find 

P(q , b 2 , b 3 , ...) = \b 2 \P( ai ,a 2 , as, . . .) • (41) 

The probability distribution of b 2 , the curvature of the 
potential U(q) at a turning point, is hence 



P{b 2 ) = \b 2 \ j ' dq db 3 db 4 . . .P(ax,a 2 ,a 3 , 



(42) 



Under the assumption that P{a\, a 2 , . . .) is free of zeros 
and divergences, this integral remains finite as b 2 — ► 0, 
and so for small b 2 we have P(b 2 ) cx b 2 . 

Small amplitude oscillations about qo have a frequency 
u) oc y/b~2, and it then follows that the probability distri- 
bution of these oscillation frequencies varies for small uj 
as 



p(u) oc oj 3 



(43) 



It is a further restriction to demand that a minimum 
at g is the global minimum of U(q). A full treatment of 
this constraint would be difficult but is fortunately not 
necessary: the crucial condition is that U (q) should have 
no nearby minima deeper than the one at qo- For that it 
is sufficient to truncate the expansion 



U(q) = U(qo) + b f(q-qo) 2 
+ fi(9-<Zo) 



9o J 



(44) 



at 0(q — qo) . After truncation, we require 1 2>3 j < \/3b 2 b4 : 
in order that U(qo) < U{q) for all q. This further sup- 
presses the probability density for small curvatures, giv- 
ing 



dqodb^... db 3 P(ai,a 2l . . .) (X b 2 . 

In turn, this brings the frequency distribution at small uj 
to the form 



p(u)) oc uj 4 



(45) 



One expects that the result of including higher order 
terms in Eq. H44J1 , and of ensuring that no more distant 
minima are lower than the one at qo, will be to change the 
constant of proportionality but not the power in Eq. Q45fl. 

Clearly, a serious limitation of this discussion it that it 
is limited to a system with a single coordinate and conju- 
gate momentum. As with our closely-related discussion 
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of random matrices, preceding Eq. (|22|l . we expect the re- 
sult to apply quite generally, provided excitations are lo- 
calized with a finite localization length at low frequency. 
In these circumstances, the coordinate q is interpreted as 
being the relevant degree of freedom for a low-frequency 
excitation within a localization volume. One of our objec- 
tives in Sec.lVland Sec. lVIl is to provide detailed evidence 
for the more general relevance of Eq. 14511 . 



V. DISCRETE ID RANDOM FIELD XY MODEL 



In this section and Sec I VII we study the excitations of 
the one-dimensional, classical XY spin chain in a ran- 
dom field. This model provides a simple but non-trivial 
example of a system with bosonic excitations which are 
not Goldstone modes, and has been discussed previously 
in Refs. and In what follows we apply the gen- 
eral approach described in Sec. II I II First, we set out 
definitions and write down the Hamiltonian TL for small 
amplitude excitations. Second, we find a local Q which 
satisfies 7i — Q T Q. As a result, we map our problem 
onto a well-known one, involving a one-dimensional chi- 
ral Hamiltonian. We next review established results for 
such one-dimensional chiral Hamiltonians, which can be 
in one of two regimes, depending on the details of their 
disorder distribution. We use the mapping to obtain p(oj) 
both numerically, and (in Sec. IVIfl analytically. 



Definitions 



The Hamiltonian for the random field XY spin chain 



is 



H 



1 N 

-Y 

2 ^ 

i=l 



uf - ^2 cos {4>i - fa+i) - Y h% cos (^ ~ Xi) 

i i 

(46) 

Here Hj are momenta conjugate to the spin angles fa, 
exchange energy is represented by — cos(fa — fa+i), and 
hi and Xi are the- amplitude and phase of a random field. 
It is convenient to introduce the notation I (fa) = — cos(fa) 
and hi(fa) = -hi cos(</> - Xi)- 

While the kinetic energy is quadratic in the momenta 
H;, the potential energy is strongly anharmonic in the 
coordinates fa. We want to find the ground state spin 
configuration fa( and the frequencies of oscillations about 
that ground state. The ground state configuration satis- 
fies 



OH 

dfa 



0. 



(47) 



Expanding H about fa to quadratic order, Eq. ll4^|) 
reduces to an expression of the general form given in 
Eq. Q, and specified by the matrices C, M and K. In 
this case, C = and M = 1. The symmetric matrix of 



spring constants, = d 2 H/dfa 
with no-zero entries 



y J 1 4>=4>o 



is tridiagonal, 



and 



dpi 



fa°_i) + dll(fa° - fa° +1 ) + dlhi(fa) (48) 



K 



= K 



= -dlKfa - 



(49) 



As discussed in Sec. IIIII with M, C and K of this form, 
the excitation frequencies of the system are Wj = 
where the stiffnesses the eigenvalues of K. Our 

tasks, then, are the linked ones of determining fa? and 
diagonalizing K. 



B. Mapping onto a chiral problem 

Since K+j is a real positive matrix, it can be writ- 
ten as the square of another real matrix Q, in the form 
K = Q T Q. Our strategy is to find Q and then study the 
related chiral Hamiltonian 



H 







o 



(50) 



As the frequencies coincide with the square roots of the 
eigenvalues of K, we will use ui to denote the eigenvalues 
of H in this and the following section. 

Because K is tridiagonal, the matrix Q can be cho- 
sen bidiagonal, with non-zero elements Qn = Ai and 
Qa-x = —Bi—i which satisfy 



Ku = Af + Bf , Ki.i+t 



-Ai + \Bi 



(51) 



To show this and find Ai and Bi, it is helpful to use the 
idea of a partial energy, first introduced in this context 
by Feigelmani^!! £i(fa) is defined to be the ground state 
energy of a subsystem which consists of sites j < i, con- 
sidered as a function of the orientation of the boundary 
spin, fa. Formally 



£ l (fa) = min V" [I (fa 
<t>3,3< 1 ^rr. 



»3+V 



i)] + hi 



(52) 



It satisfies the recursion relation 



£i(fa) = min [/ (<j> - i>) + £_i(V0] + hi(fa) . (53) 

Now let the value of tp which results from the minimiza- 
tion in Eq. (|53|l be rpo(fa), and define fa—\ as a function 
of fa by fa—i(fa) = ipo(fa)- With this notation, the con- 
dition that the right side of Eq. i|53|) is at a minimum 
takes the form 



■d A I{ 



fa-i) + d$£. 



i-i W-i 



)=0. 



(54) 



By differentiating Eq. I|54|) with respect to fa , remember- 
ing that fa— i is a function of fa in the sense described 
above, we find 



dll(fa-fa-!) 



din 



fa- l ) + d 2 E l -i(fa- l ) 



(55) 
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In addition, we differentiate twice with respect to 
<pi, again remembering that (f>i-i is a function of 4>i, to 
find 



i-l 



(56) 



This allows us to solve Eqns. ijIFfy. (gHJ and (f5*T)) . obtain- 
ing 

2 



At 



din 



din 



die, 



and 



5/ = 9iJ( 



-$) + afc($) : 



(57) 



(58) 



which completes the derivation of Q. 

Let us examine the chiral Hamiltonian Ti, Eq. (|50() . By 
rearranging its rows and columns it can be put into the 
form 



n = 



( 


Si 








. 


. 


\ 


Bl 





^4 2 





. 


. 








A 2 





B 2 


. 


. 











B 2 





^3 • 


. 














^3 


. 


. 

















. 


. 


Bn 


V o 











. 




o 1 



(59) 



which is familiar as a one-dimensional tight binding 
model with only off-diagonal disorder and is also 
referred to as a random chiral one-dimensional Hamilto- 



and the localization length of these states diverges for 
uj — ► as 



£(w) = |logH| 



(62) 



By contrast, for a / (sometimes referred to as the 
staggered regime^) the density of states varies as 23 



p{uo) oc uo 13 , uj ^ a, 



(63) 



with a power /3 that depends on a and a. The localization 
length is finite in the small u> limit, and varies for small 
a as £ oc \a\ . 

While the Ai and Bi which arise in our treatment of 
the XY model have specific correlations not present in 
the chiral problems studied previously, some comparisons 
are useful. In particular, considering for simplicity weak 
disorder, hi <C 1, the system arising from the XY model 
turns out to be in the staggered regime. To see this, 
note that for weak disorder, 

dH(4>9-d>9,)^l, 



< 1, so that 



Ai a 1 - ^Si-Mti), B< « 1 + (64) 

and hence a = (d^£)/2. We argue in Sec. IVII that 
> 0. We therefore expect low frequency states 
to be localized with localization length £ oc and 
with a power law density, as in Eq. I|63|l . Quite sepa- 
rately, if states are localized, the assumptions that led to 
result p{lo) oc a; 4 in Sec. II VI are justified, and so we ex- 
pect the exponent (3 = 4. We postpone further analytical 
work to Sec. IVII and first treat the problem numerically. 



C. Established results for one-dimensional chiral 
problems and implications for XY model 

There has been extensive previous work on one dimen- 
sional models of the type represented by Eq l|59(l . with 
disorder in the Ai and Bi which is uncorrelated and cho- 
sen to have a simple, known distribution. The results 
serve as a guide for our calculations, and we summarise 
them here. Parameterise the matrix elements Ai and Bi 



Ai = 1 



SAi, Bi = l- 



5Bi 



(60) 



and take 5Ai and SBi to be independent Gaussian ran- 
dom variables with zero mean and standard deviation a. 
Behavior at small lu is very different according to whether 
or not a is zero. 

For a = and u <C o, the density of states has a 
singularitji^i at ui = of a type first obtained in a related 
problem by Dyson, 



p(uj) oc 



1 

to I log 3 uj I 



(61) 



D. Numerical study of the one-dimensional 
random field XY model 

Our numerical procedure involves several steps. First, 
for a system of length L, we construct the function £i(4>) 
for each i by iterating Eq. (|53|l from i = 1 to i = L. 
Then we determine (jjf , iterating from i = L to i = 1 and 
using the fact that for each 0°, ffi-i minimises the right- 
hand side of Eq. (|53|) . Knowing the ground state spin 
configuration, we compute the matrix elements appearing 
in the chiral Hamiltonian Eq. I|59|l . using Eqns. (|57|l and 
(|58|l . Finally, we employ the transfer matrix technique 
developed specifically for such Hamiltonians in Ref . l24l to 
find the integrated density of states. For the random field 
we choose a uniform distribution of [hi cos(xi), hi sin(xi)] 
over a disc of radius D, independently for each i. We find 
that it is sufficient for each disorder strength D to study 
a single realization in a system of length L = 10 6 . 

Our results for the integrated density of stiffnesses, 
N(ui) — J Q duj'p(uj'), are shown in Fig 1. Behavior at 
large uj approaches that in the disorder-free system, as 
indicated. At small uj we expect the power law 



N(w) oc uj {fM 



(65) 



10 



10° 



produce a universal density of states, p(u>) oc u 4 . In the 
next section we present the analytic derivation of this 
result. 



3 



10' 



10' 



10" 



10" 



10" 



10 u 



FIG. 1: The integrated density N(lo) plotted as a function 
of frequency lo using logarithmic scales. Dashed, dotted and 
dot-dashed lines represent disorder strengths D — 0.3, 0.1 
and 0.01, respectively. For all three cases, the integrated den- 
sity converges at large ui to that of the disorder-free system, 
represented by the full line which has gradient 1. 



This form and also the value p = 4 are both already ap- 
parent in Fig 1, as was also the case, at lower precision, 
in earlier numerical results of FogleuSi for much shorter 
chains. It is more instructive, however, to observe that 
for a power-law density, dimensional analysis fixes the de- 
pendence on disorder strength D to be p(u>) cx D~ 2 ^^uj^ . 
In Fig. 2 we plot D 8 / 3 N(oj) as a function of w, showing 
collapse of data at small oj for three different disorder 
strengths, and power-law behavior with [3 — 4. 



10° 



3 

£ 10° 



10" 



10" 



10" 

CO 



10 u 



FIG. 2: D a/3 N(oj) vs oj, for disorder strength D = 0.3, 0.1 
and 0.01, shown with dashed, dotted and dash-dotted lines, 
respectively. The straight line has a slope (/3 + 1) = 5. 

We conclude that, while the exponent f3 appearing in 
the density of states is, for a generic chiral problem, 
disorder-dependent and hence non-universal, the partic- 
ular disorder generated in the mapping from the ran- 
dom spin chain, following Eqns. (|53|l . I|57|) and l|58|) . has 
a specific distribution and the correlations necessary to 



VI. ANALYTIC TREATMENT OF THE 
RANDOM FIELD XY MODEL IN THE 
CONTINUUM LIMIT 

In this section we continue our examination of exci- 
tations in the one-dimensional random field XY model, 
using the continuum limit to make analytic progress. We 
find both the localization length as a function of disor- 
der strength and the density of states as a function of 
frequency, confirming heuristic arguments and numerical 
results given above. 



A. Mapping to a chiral problem 

The continuum limit of the random field XY spin chain 
is reached at weak disorder, hi <C 1. In this limit it is 
possible to replace the discrete index i with a continuous 
variable x, and the Hamiltonian of Eq. I|46(l becomes 



H 



dx 



1 



-U 2 



(d x 4>) + h((j>(x),x) 



(66) 



with h((p(x), x) = —h(x) cos (</>(x) — x( x ))- As in the dis- 
crete case, we are interested in the configuration of spins 
4>o(x) that minimizes the potential energy. This configu- 
ration satisfies 



dl<j> (x) +d <i> h(<f>o (x),x) = 



(67) 



and the amplitudes of normal mode excitations about the 
ground-state 4>o obey 



-dl + dlh{<t> {x),x)]^ = w^. 



(68) 



Just as in the discrete case, because this equation de- 
scribes deviations from a minimum, all normal modes 



have positive stiffness k 



As a result, the operator 



appearing in Eq. (|68J) can be written as a square. We set 
d 2 



dx 2 



+ d<t>H<Po(x),x) 



Q Q, 



(69) 



where Q = —d/dx + V(x). In other words, we require a 
function V(x), which we term the chiral potential, that 
satisfies 



dV(x) 
dx 



+ v 2 ( x ) = d 2 h(Mx),x)- 



(70) 



In order to understand properties of the chiral poten- 
tial, it is useful to introduce £{<fi, x), a continuum version 
of the partial energy^ defined following Eq. I|52|) by 



£(ip,x) 



mm 

b(x)=<P Jo 



dy[{d y <t>) 2 + h{4>{y),y)\ • (71) 
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Similarly, the continuum version of Eq. (|53|) is 



d x £ 



h((j>,x) . 



(72) 



which can be thought of as a Hamilton- Jacobi equation 
for the action £ of a particle with coordinate <f> moving 
as a function of time x in the time-dependent potential 
h(4>(x), x). In addition, the continuum version of Eq. I|54|) 
relates the ground state configuration to this action, via 



d(f>o (x) 
dx 



d<j,£((f>Q(x),x) 



(73) 



(74) 



It is easy to check that by writing 

V{x) = dl£{4>v{x) : x) 

we solve Eq. {7UJ). Thus the chiral potential V(x) may 
be expressed simply in terms of the dependence of the 
ground state energy of a half-system on the boundary 
spin, 4>{x): the half-system has coordinate y taking values 
in the range < y < x. With V(x) in hand, we wish 
to study the continuum version of Eq. 159|) : the chiral 
operator 



H = 



Q 
Q T 



dx 



o 

VV{x) 







(75) 

which has eigenvalues ±a-> that are the square roots of 
those appearing in Eq. JHSJ). 



that, by assumption, V(x) is mainly positive, and for 
positive V(x) Eq. l(77|) has stable fixed points close to 
9 = rwr, with n integer. Rare fluctuations of V(x) which 
are negative for a long interval in x allow 9(x) to grow, 
evolving with increasing x from one such fixed point to 
the next. Suppose 0(x) leaves the vicinity of one fixed 
point at x = x\ and arrives in the vicinity of the next at 
x — X2- For xi < x < X2 we neglect u in Eq. I|77|l and 
obtain 



(W 



sin (26») 



dx V(x) . 



(78) 



We estimate the integral on the left hand side by noting 
that the most important contribution comes from the 
end points, where Q(xi) n7r, 9(x 2 ) ~ (n + l)n and 
sin (29) ~ u). For 9 to increase by tt we therefore require 
a fluctuation in V(x) which is sufficiently negative and 
extends over a sufficiently large interval in x that 



log (uj) > dx V(x) . 



(79) 



Let P(u>) be the probability per unit length for such a 
fluctuation to occur: the integrated density of states, 
Jq p(ui)duj, is simply P(t-o). It is natural to expect this 
probability to be exponentially small in log(w) for small 
lu, so that, introducing a constant a, 



P{oj) ~ exp(a log(w)) = uj c 



(80) 



B. Treatment of one-dimensional chiral problems 

As with the lattice version, discussed in Sec IV CI the 
spectral properties of H, Eq. i(75|) . have been studied ex- 
tensively with simple choices for the probability distri- 
bution of V(x). Behavior is as summarised for the lat- 
tice version in Sec IV CI a particularly detailed study is 
given in Ref. El Here, for completeness we sketch the 
derivation of the result that is of most importance for 
our work: the density of states at low frequency in the 
staggered regime, where (V(x)} > for the continuum 
system plays the same role as a > for the lattice model. 

Following Ref. l4ll consider the coupled first-order dif- 
ferential equations 



and hence 







V(x) 







ipi(x) 
ip 2 (x) 



ipi(x) 

V>2(z) 



(76) 

for x > 0, with boundary conditions ifti(0) = 0, ^2(0) = 
1. From the node counting theorem, the integrated 
density of states is given by the density of zeros of 
ipi (x) per unit length. Introducing the parametrization 
-01 {x) — p(x) su\9(x) 1 '4>2(x) — p(x) cos9(x), one has 



d9(x) 
dx 



= ui- V{x) sin (20(x)) 



(77) 



Thus we require the average rate of increase in phase 
with length, d9(x)/dx. To find this at small oj, note 



p(uj) 



(81) 



with (3 — a — 1. From an extension of this approach, one 
also finds^i that the localisation length £ at low frequency 
varies with the staggering (V(x)) as £ ~ (^(.t)) -1 . 



C. Calculation of the density of states 

A central difficulty of our problem, of course, is that 
V(x) does not have a simple, given distribution. Instead, 
it should be determined by solving Eq. (|70|l . or calcu- 
lated from the partial energy using Eq. I|74|) . after in 
either case first finding the ground state configuration 
(j)o{x), using Eq. (|S7Jl or Eq. J73J). Our detailed calcula- 
tions are based on Eqns. H7U|) and Ht)7fl . Before presenting 
these calculations, it is useful to develop some qualitative 
understanding by an alternative route, using the partial 
energy, £(x,4>), and its connection to Burgers turbulence. 

The evolution of £(x,4>) with x is described by 
Eq. I|72|) . which is similar in form to the much-studied 
KPZ equation^ In this correspondence, x and <p play 
the roles of time and space coordinates, respectively. The 
stochastic KPZ equation, however, which reads 



d x £ + \ (d^£f 



Ddl£ = h{<t>,x), 



(82) 
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differs from Eq. (|72|l in two respects. One is in the ab- 
sence of the dissipation term, with coefficient D; the 
other is in the correlations of the force, h((/),x). The ab- 
sence of dissipation is of limited importance, because the 
relationship between the solutions of Eq. I|82[l in the small 
D limit and those of Eq. (|72(l is well-understood: while 
Eq. (|72|) generally has solutions with many branches, cor- 
responding to local minima in the energy of the spin 
chain, by taking the solution of Eq. (|82|l for D — > one 
finds the envelope of absolute ground states as a function 
of the boundary spin <fi(x) . For this reason, it is rather 
natural to study Eq. (|S2*|l in our context. 

By contrast, the nature of force correlations is more 
significant: while in standard form the KPZ equation has 
a random force h((f>, x) that is white noise in both x and 
0, our interest lies with correlations that are long-ranged 
in 4>, and have the form 

( h{<t>i,Xi)h{<t>2, X2) ) = hcos(4>i - <fa)6(xi — x?) . (83) 

Such correlations have been studied previously, in the 
context of Burgers turbulence, where the Burgers equa- 
tion, Eq. (|HIJ), (equivalent to Eq. with D — ► 0) to- 
gether with the correlator Eq. (|83|l . describes motion of 
a one-dimensional fluid. 

In this analogy, we view 4>(x) as the coordinate of a 
particle as a function of time, x, and interpret Eq. i|t)7|) 
as the equation of motion for the particle. Imagine a 
fluid of such particles, without pressure and stirred ran- 
domly with force correlations derived from Eq. (|83|) . Let 
w(</>, x) = d(j)/dx be the velocity, which satisfies the Burg- 
ers equation 

d x u + ud^u = d^h{(t),x) , (84) 

where comparison with Eq. I|72|l shows that u — d^S. 
The chiral potential is therefore given by V(x) = 
d < pu((f>o(x), x), the velocity gradient in a one-dimensional 
fluid flowing in a time-dependent potential h(cj),x), cal- 
culated at a point that moves with the fluid. 

From literature on the Burgers equation (see Refs.l43l 
I44TI45F) . or alternatively by thinking about the ground 
state of a spin chain as the boundary spin 4>(x) is varied^ 
one arrives at the following picture for £((/), x). As a func- 
tion of <p it typically consists of a few piecewise smooth 
sections, which meet at cusps that are local maxima of 
£{4>,x). These cusps in £(cj),x) are negative discontinu- 
ities or shocks of the Burgers velocity field, u{4>,x), at 
which 

lim \u{4> + e, x) - uU, x)] < . (85) 

They occur at the points where the ground-state spin 
configuration of the half-chain changes discontinuously as 
</> varies. With increasing system length x, they undergo 
an evolution in which existing cusps merge and new cusps 
are born, at matching average rates. If the trajectory of 
a particle moving with the Burgers fluid should meet a 
shock, the particle remains trapped by the shock for all 



subsequent x. From these statements it is clear that (for 
almost all boundary conditions 4>(L)) the ground state 
trajectory <fio(x) does not intersect any shocks. As a next 
step, from this we expect that (d^£(4>o(x), x)} > 0, on the 
basis that, first, an average of d^£((j), x) = d^u^, x) over 
all 4> must be zero, since u(<f>,x) is periodic in <j>, while, 
second, an average restricted to (f>o avoids negative dis- 
continuities of it(<^>, x). Finally, returning to excitations 
of the spin chain viewed as a chiral problem, we con- 
clude that this is in the staggered regime and expect a 
finite localization length £ ~ (d^£((f)Q(x),x)) > 0; from 

dimensional analysis, we expect £ ~ h^ 1 ^ 3 . 

Two weaknesses of the argument we have given are 
clear: first, a more detailed treatment of averages over 
ground states (f)o{x) would be desirable; and second, it is 
not certain that behavior known for chiral problems with 
disorder uncorrelated in x will necessarily be present in 
our system, with correlations of V(x) determined from 
ground-state properties. Nevertheless, results of the de- 
tailed calculations below are consistent with the foregoing 
conclusions. 

To make further progress, we need to study statistical 
properties of V(x). We find that a direct attack on this 
problem using the boundary conditions that are physi- 
cally appropriate (in which the values of <po{ x ) are spec- 
ified at x — and x — L) is too difficult mathematically. 
Instead, we approach it indirectly, by relating it to a ver- 
sion of the problem with simpler boundary conditions, in 
which 4>o(x) and d x (j>o(x) are specified at x = 0. With the 
latter boundary conditions, we solve jointly Eq. (|67|) for 
4>o(x) and Eq. i|7U)l for V(x). To discuss the connection 
between the two alternative sets of boundary conditions, 
consider Eq. (|70|l with an arbitrary forcing term, written 
as f(x): 

*^- + V*(x) = f(x). (86) 

This equation should be integrated from x — towards 
x = L, with initial condition V{x) — > oo for x — > since 
this is the behavior of d^£((j)o(x),x) at small x and V(x) 
is related to the partial energy by Eq. J73J. Crucially, 
Eq. H86[) is unstable in the sense that, if V(x) reaches 
sufficiently large negative values for V 2 (x) to dominate 
over f(x), the solution escapes to towards V{x) = — oo, 
with the form V(x) cx (x — xq)" 1 as x approaches Xq, 
the position of the instability, from below. If such an 
instability is reached, it signals the fact that the tra- 
jectory <f>o (x) no longer represents the absolute ground 
state of the spin chain but is instead a local maximum 
in the energy as a functional of configuration. If the 
forcing term is derived from the ground state configura- 
tion using f(x) — dj ) h(^)o(x),x), it will have correlations 
that ensure this instability is never reached. But if we 
treat Eq. H86fl a s a Markov process in the way we set 
out below, some realizations will prove unstable in this 
sense. Such trajectories should be discarded, and this 
can be arranged by supplementing Eq. (|86|) with absorb- 
ing boundary conditions at V(x) = —oo. The surviving 
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trajectories must be weighted in order to sample ground 
states appropriately, and we return to this aspect in due 
course. 

First we check that, specifying 4>o(x) and d x <po(x) at 
x = and integrating Eqns. l|H7l) and ijTUI) together, we 
indeed have a Markov process. This is seen most clearly 
by returning to the discrete equations, taking Eq. (|56J) in 
place of Eq. (|7U)l and setting d\l w 1 to obtain 



Vi 



Vi-i 



1 + Vi. 



dlhi 



(87) 



Similarly, in place of Eq. l|73|l . Eq. I|54|) can be written as 



4>i = <f>i-l + 9<j>£i-l(<l>i-l) ■ 



(88) 



It is now clear that, since 4>i is determined by hj for 
j < i, it is independent of the function hi(4>). Moreover, 
by construction, each hj(<f>) is an independent random 
function. Returning to the continuum limit and noting 
the correlator for h(<f>, x) given in Eq. (|83[1 . we see that in 
Eq. (|51)|) we should take f(x) Gaussian distributed, with 
zero mean and correlator 



(f(x)f(y))=hS(x-y). 



(89) 



In light of the discussion given in Sec. IVI Bl our next 
step is to find the probability of an unusually large neg- 
ative fluctuation in V(x), integrated over an interval of 
length I = x% — X\. To this end, it is helpful first to define 
S(a) by the equation 



exp ^— a J dx V(x)^j ^ = exp 



(-lS(a)) (90) 



where the angular brackets indicate averaging over V, 
and we will find that S(a) is independent of I when I is 
large. We denote this probability (which appeared previ- 
ously in Eq. JHDJl) by Pq(u>) below, and calculate it using 



P (w) = max / da exp (-lS(a) + alog(w)) . (91) 



At small to we can approximate the integral by its value 
at the saddle-point, determined from 



^ + log(w)=0. 



(92) 



In turn, the maximization on I (remembering that the 
value of a at the saddle-point is itself a function of I) 
gives 



Therefore 



5(a) = 0. 



(93) 



(94) 



where ao is the solution to the Eq. I)93[l. 

Now we must compute S(a). A convenient method is 
to derive from the Langevin equation Eq. (|86|l a Fokker- 
Planck equation, make a similarity transformation of the 
latter to obtain a Schrodinger equation, and express this 
as a path integral. The absorbing boundary condition 
at V = — oo is built in automatically in this approach. 
These methods are described, for example in Ref. In 
this way we find 



exp (—IS (a)) 



J VV exp 




- V(l-a)} 


/ Wexp 


-jo** U((^) 2 + 









(95) 



The path integrals in this expression are propagators in 
imaginary time x for a particle moving with coordinate 
V (x) in a polynomial potential, which is V A — V(l — a) in 
the case of the numerator and V A — V in the case of the 
denominator. For large I, both path integrals are dom- 
inated by ground-state contributions, justifying the l- 
dependence displayed in Eq. I|90|l . In this limit, moreover, 
S(a) is given by the difference in ground-state energies 
for the two potentials. It is clear that S(a) = 5(2 — a) 
and that 5(0) = 0. Hence = 2 and 

P (w) oc w 2 . (96) 

It is incidentally also apparent that (V) > 0, confirming 
our earlier argument that the chiral description of spin- 
chain excitations is in the staggered regime. 



To complete our calculation of the density of excita- 
tions in frequency, a further step is necessary. We have 
so far considered spin configurations 4>q(x) that are gen- 
erated in the ensemble of disorder realizations using fixed 
values for (f>o(x) and d x (f>o(x) at x — 0, and are local min- 
ima of the energy but not necessarily the absolute mini- 
mum. We should weight these configurations by a factor 
Pi(w), according to the probability that they appear in 
an ensemble with physically appropriate boundary con- 
ditions, in which 4>q(x) is fixed at x — and x = L. 
In addition, to obtain properties of excitations about the 
ground state, we require a further weighting factor P2(uj), 
involving the probability that a configuration is the ab- 
solute minimum in energy. This is in principle a diffi- 
cult quantity to determine, 46 because it involves global 
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features of the system, and we content ourselves with a 
heuristic approach which is in a similar spirit to our dis- 
cussion of an anharmonic oscillator in Sec. IIV Bl Specif- 
ically, we discard all configurations that have a nearby 
maximum of the energy. In this way we correctly exclude 
local minima that are separated from the ground state by 
the nearby maximum, but we make errors of two kinds. 
First, we fail to exclude local minima that have no nearby 
maximum and are separated by a large distance in con- 
figuration space from the true ground state: we assume 
(as in Sec. IIVBII that excitations about such local min- 
ima have the same statistical properties as those about 
the true ground state. Second, we wrongly exclude true 
ground states that have nearby maxima: we assume that 
those ground states which remain are characteristic of 
the whole set. We remark finally that if the factor P^ioS) 
is omitted, we obtain an excitation density averaged over 
all configurations that are local energy minima. 

To find the (^-dependence of these two weighting 
factors, we consider a family of nearby configurations 
4>{s,x) = (f>o(x) + r)(s,x) for x\ < x < x 2 , parameter- 
ized by s. The corresponding family of chiral potentials 
is V(s, x) — V(x) + W(s, x), and we restrict attention to 
small r)(s, x) and W(s, x). The weight -Pi(w) appears be- 
cause, in disorder realizations which generate negative 



fluctuations of V(x), trajectories of <j)(s,x) as a func- 
tion of x are compressed by an amount which increases 
with the size of the integrated potential fluctuation. It is 
therefore determined by comparing rj(s, x%) with rj(s, x\). 
The weight Pa(^) is determined by finding the probabil- 
ity for escape of W(s,x) to negative infinity, signaling 
the occurrence of an energy maximum. 

We find the evolution of rj{s, x) and W{s, x) in terms 
of V(x) by linearizing Eqns. i|t)7[l and (|7L)|) respectively, 
obtaining 



drj 

— = V{x)r]{s,x) 



(97) 



and 



dW 



— + 2V(x)W(s, x) = d$h(4>, x)r)(s, x) . 



These equations have the solutions 



rj(s, x) — exp 



dy V(y) 



r)(s,xi) 



(99) 



and 



W(a,x)=exp[-2 dy V(y) \W{s, Xl ) + v{s, Xl ) dy 



exp 3 / dzv(z) d^h(My),y) 



(100) 



Let us choose the parametrization s in such a way that 
at s = 0, r)(s,Xi) — and V(s,Xi) has a minimum. 
Then for small s, rj(s,xi) oc s and W(s,xi) oc s 2 . We 
find P\(oj) as follows. For the integrated value of V(s, x) 
along a trajectory from the family to have a value similar 
to that at s = 0, the integral of W(s,X2) must not be 
too large. We hence select those values of s for which 
f x l dxW(s,x) < 1. Using Eq. ljTUU|) and that fact that 

J xi 2 dxV(0, x) = log(w), this requires |s| < cj, which in 

turn implies that |^(s, xi)\ < u> and \r)(s, x 2 )\ ^ w 2 . Since 
the physically relevant boundary conditions for the spin 
chain fix the value of (j)( x ) &t x = L (that is, at larger 
x), we weight configurations uniformly in rj(s,X2) and 
conclude that -Pi(w) oc uo 2 . Turning to ^(w), we require 
that W(s,X2) calculated using our linearization should 
not be large and negative for any s, since otherwise a 
full treatment, including non-linearities, would with high 
probability result in escape. Minimizing W{s,X2) with 
respect to s, using Eq. (|100f> . we require 



min W(s, x 2 ) <x 



-lo- 2 I 2 



> 



1. 



(101) 



where 



d:X 



exp 3 dyV(y))d$h(Mx),x) 



(102) 

For small ui and large \x 2 — X\\, / is a random variable 
whose distribution is independent of w, and hence the 
probability that min s W(s,X2) ^ —1 is P2(oj) oc ui. 

Combining these results with the form derived for 
Pq(uj), we find an integrated density of states for excita- 
tions, averaging over all local minima of Eq. (|46|l . which 
varies as 



N(w) oc P (w)Pi(u) = lu 4 



(103) 



while for the global minimum we find 

N(lo) oc P (lu)P 1 (lu)P 2 {uj) = uj 5 . (104) 

Hence = 3 or /3 = 4, in agreement with Eqs. ll4^|) 
and l|45() . Strikingly, the behavior derived here for the 
random field XY chain matches that expected from the 
simple discussion of an anharmonic oscillator, given in 
Sec lrVBl 
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VII. OTHER SYSTEMS WITHOUT 
GOLDSTONE MODES 



and may be used to find the ground state configuration 
via the analog of Eq. J73J), 



In this section we discuss the extent to which the ap- 
proach we have set out for the random field XY spin 
chain can be extended to treat excitations in other mod- 
els for disordered systems without a broken continuous 
symmetry. The XY spin chain has two obvious and im- 
portant simplifying features: disorder couples only to one 
of the dynamical variables (<f> but not II) , so that frequen- 
cies are simply related to stiffnesses; and the system is 
one-dimensional. The calculations we have described in- 
volved several steps: mapping to a chiral formulation; 
analytic determination of the excitation density in fre- 
quency via a study of statistical properties of disorder 
within this chiral formulation; and, for an efficient nu- 
merical treatment, the use of a partial energy. As we 
broaden the range of models under consideration, fewer of 
these steps remain possible. We consider, first, a generic 
one-dimensional continuum system, with a pair of conju- 
gate dynamical variables at each point and disorder that 
couples to both, and second, a specific example of such 
a problem, the random field Hcisenberg spin chain. In 
these cases we find a chiral description and introduce a 
partial energy, but are not able to determine statistical 
properties of the disorder appearing in this description. 
Instead, we derive for such problems the general relation 
between densities of stiffnesses and of excitation frequen- 
cies, suggested from an analysis of random matrix theory 
in Eq. I122H above. Third, for the random field XY model 
in two dimensions we show how to introduce a chiral de- 
scription, but leave applications of this for future work. 
For all these problems, we expect p(u)) (x w 4 at small u>, 
on the general grounds discussed in Sec. IIIII 



A. Generic one-dimensional problem 

Consider a one-dimensional system with conjugate dy- 
namical variables 4>i(x) and cj>2(x), and the Hamiltonian 



H = dx 




(105) 

We wish to study harmonic excitations about the ground 
state: for i = 1,2, let </>° be the configuration that min- 
imises the energy, Eq. I|105|) . and write = $ + tpi. 
Expanding to quadratic order in 

-9x + 9lh dfadfrh 



H 



(106) 



To write this as Ti — Q T Q and find Q, we essentially 
repeat the sequence of arguments which led from Eq. I|66|) 
to Eq. (175)1 . We introduce a partial energy £(<fii,(p2,x) 
which satisfies 



d x £ 




Kfafax) (107) 



(108) 



Then we define a 2 x 2 matrix version of the chiral po- 
tential, via 



(109) 



It satisfies 



dx 



Vij + V ^ V kj = dtod^hifafax) (110) 



fc=l 



and therefore Q can be chosen to have elements 



V, 



ij ■ 



(111) 



In this way, we have written Eq. (jl()6J) as a type of chiral 
problem, which we term two-channel because V is a 2 x 2 
matrix. 

Further analysis can be divided into two stages. One 
stage, given V, is to find d{n) and p(u>), as was done 
for the single-channel problem in Sec. IVIBI The other 
stage is to determine the distribution for V, as was done 
for the single-channel problem in Sec. IVI CI General 
one-dimensional multichannel chiral problems of the type 
that arise in our calculation of stiffnesses have been stud- 
ied in Refs. |2(| and 0, with V chosen Gaussian dis- 
tributed and uncorrelated in x. However, as far as we are 
aware, there has been no previous work on multichannel 
chiral problems of the type that give frequencies. For 
the two-channel problem, while we have not been able 
to make progress in obtaining the distribution of V, we 
have been able to find a general relation between calcu- 
lations of stiffnesses and frequencies. This connects 
and p(lo) in the way given by Eq. I|22[) . We present these 
arguments next. 



B. Two-channel Bosonic Chiral Problems 

In this subsection we compare the calculation of stiff- 
nesses k = A 2 , determined from the eigenvalue problem 











U<2 



x 



(112) 



with the calculation of frequencies uj = e 2 , determined 
from 











W2 J \0-2W2 



(113) 



We do this following the technique set out for the single- 
channel problem after E q. 117611 . and described for multi- 
channel problems in Ref.|48l For both stiffnesses and fre- 
quencies, we write ipi = ai/ , 2, where a is an x-dependent 
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2x2 matrix. The eigenvalue density, d(n) or p(u>), is 
then determined from the evolution of a with x, via a 
node-counting theorem. The evolution equation for a is, 
in the case of stiffnesses, from Eq. Q112J) . 



dx 



+ a 2 s ) - a s V - V T a s 



and in the case of frequencies, from Eq. (|113|> . 



daj_ 

dx 



(a 2 + a 2 f ) - a f V -V T a } . 



(114) 



(115) 



These equations are the equivalent for the two-channel 
problem of Eq. I|77(l for the single-channel problem. Fol- 
lowing Ref. [4a the density of states is given by the rate 
at which the eigenvalues of a s or a/ move with increasing 
x along the real axis. Our aim is to compare this rate in 
the two cases, taking k and uj small, in a way that does 
not require detailed knowledge of the distribution of V. 
We assume only that (as for single-channel problem) V 
fluctuates about a non-zero mean, and that both the fluc- 
tuations and the mean of V have comparable importance 
in the evolution of a s and a/ with x. 

In outline, this evolution is as follows. If fluctuations 
in V are omitted, a s and a/ have stable fixed points at 
which both their eigenvalues are small in k and u> respec- 
tively. In both cases, the stable fixed point has a basin 
of attraction with a boundary that is reached when an 
eigenvalue of a s or aj is large and positive (C(A _1 ) or 
0(e~ 1 ) respectively). Restoring fluctuations in V, with 
increasing x we find (partly on the basis of simulations, 
not presented here) that a s or a/ fluctuates in the vicin- 
ity of its fixed point for intervals that are long if n or 
lj are small. Each such interval ends when fluctuations 
take the matrix to the boundary of the basin of attrac- 
tion. One eigenvalue of o s or o/ then runs off to positive 
infinity, reappears from negative infinity and returns to 
the vicinity of the fixed point. This is entirely analogous 
to the evolution in the single-channel problem of tan(0) 
as 9 increases from 9 w nil to 9 w (n + l)7r, which is 
described in Sec. I VI Bl 

To develop the picture further, we consider separately 
the regions close to the fixed points, which are different 
for a s and af, and the region far from the fixed points, 
which is essentially the same in both cases. It is useful 
to introduce an explicit coordinate system. We write 



^1 + 6^1 +603): 



(116) 



where £0, £1 and £3 are random with mean zero, and (V) 
is taken proportional to the unit matrix, with a propor- 
tionality constant that can be changed by a rescaling of 
the length coordinate x. We also set 



a s j = si + z\o\ + z 2 a 2 + Z3CT3 



(117) 



For calculation of stiffnesses using this coordinate sys- 
tem, one has z 2 = and Eq. (|114|) becomes 



dz\ 
dx 



^p- = 2Xsz 3 - z 3 + £ Q z 3 + £ 3 s , 
dx 

^- = X{s 2 + z 2 + z 2 )+X-s + i Q s + + £32:3 . 

If fluctuations in V are omitted (by setting = for 
i = 0, 1 and 3), these equations have a stable fixed point 
at s w A, z\ = z 3 = 0. Including fluctuations in V, the 
typical magnitudes of s, Z\ and z 3 are 0{\). 

By contrast, for calculation of frequencies Eq. I|115|) 
gives 

dz\ 

— = 2esz x - Zi +£ zi +€is, (119) 
dz2 

—— = 2esz 2 + e - z 2 + Ko z 2 , 
dx 

dz 3 

— = 2esz 3 - z 3 + £ a z 3 + £ 3 s , 
dx 

J- = e (s 2 + z 2 + z 2 + Z f) - s + Cos + £1*1 + 6^3 ■ 

In this case, if fluctuations in V are omitted there is a 
fixed point at s ~ e 3 , z 2 ~ e, Z\ — z 3 — 0. Including 
fluctuations, the typical magnitudes of s, z\ and z 3 are 
0(e 3 ), while z 2 remains 0(e). 

Now consider escape of a s or a/ from the vicinity of the 
relevant fixed point, which requires s X in the case of 
stiffnesses and s ^ e 3 in the case of frequencies. In these 
regimes we argue that the evolution equations in the two 
instances are essentially equivalent. More specifically, we 
make two approximations. First, we take z 2 — in both 
cases, even though this is exact only in the first case. We 
do so because in the second case it is clear from Eq. <|119f) 
that z 2 ~ e provided s < (2e) _1 , and hence that non-zero 
z 2 has no important effect on the evolution of z\ , z 3 and 
s. Second, for a s or a/ far from its fixed point, we omit 
the terms independent of a s or aj (and small in k or uj) 
from the right-hand sides of Eqns. (|114[) and l|115|l . With 
these approximations, and making the rcscalings a = Xa s 
and a — eaf, the two equations both become 



— = a - aV - V a. 
dt 



(120) 



2Xszi — z\ + ZqZx + £is . 



This stochastic process results in a stationary probabil- 
ity distribution for a if absorbing boundary conditions 
are imposed at infinity and probability flux is injected 
near a = 0. In the stationary state there is an eigen- 
value flux along the positive real axis, with a rate that 
is determined by the probability density in a region near 
a — 0. The size of this region, and therefore the prob- 
ability density within it, follows in each case from our 
discussion of the fixed points and their neighborhoods. 
Characterizing a by the value of its trace, and combining 
the fixed-point coordinates with the rescalings used to ar- 
rive at Eq. 1120f) . the value of Tra in this region is 0(X 2 ) 
in the first instance and C(e 4 ) in the second. From this, 
we conclude that the integrated densities of stiffnesses 

(118) 

and frequencies are equal for A 2 ~ e . We hence recover 



17 



from this discussion of one-dimensional systems the re- 
sult given in Eq. I|22l) . which was reached in Scc. lllll bv a 
very different route. 

For completeness, we remark that the derivation pre- 
sented here crucially depends on (V(x)) ^ 0, and as- 
sumes that £o, £i; and £3 become comparable to 1 
only during their rare large fluctuations. And indeed, 
for ( V(x) ) = it has been conjectured and checked 
numerically^ that in that case p{ui) oc w" 1 / 3 even though 
it is well knowni^ that d(k) oc log(«;). 



C. The random field Heisenberg spin chain 

The random held Heisenberg spin chain is perhaps 
the most obvious example of the generic one-dimensional 
problems discussed above. The conjugate dynamical 
variables are, of course, the two components for displace- 
ment of a spin from its orientation in the ground state, 
and phase space is made up of spheres for each point on 
the chain. The curvature of phase space is responsible 
for some changes in formuls derived in Sec. IVII Al which 
we now set out. 

The continuum version of the random held Heisenberg 
spin chain, parameterizing spin orientations using the an- 
gles 9{x) and <fr(x), is 



H 



dx 




+ sin 2 (6>) 



dx 



h{4>, 9, x) 



The evolution equation for the partial energy is 

h(cj),e,x) 



d£ 1 { fd£\ 2 1 /<9£ x r 



dx 2 \\dOJ sin 2 ((9) \d(j) 



and, given £, the ground state configuration can be cal- 
culated using 



-o&t , — — = dot 



D. Two-dimensional random field XY model 

The two-dimensional random field XY model provides 
an example for which our formulation of excitation prob- 
lems using chiral Hamiltonians can be carried through in 
more than one dimension. The two-dimensional version 
of Eq. ((SSJl is 



K = 



n 2 + {d x <f>) 2 + (d v fif + h(<t>, x, y) \ dxdy 



(121) 



(122) 



and the ground state 4>o satisfies 

" d 2 d 2 1 

" 0^2 + 0^2 ^0 +d <l> h(<f>o, x, y) =0. 

The amplitudes of normal mode excitations obey 



d 2 d 2 



We introduce a chiral potential V(x,y), defined as the 
solution to 

d 2 x V + (d x V) 2 + d 2 V + (dyV) 2 = 8lh(<j> ,x, y) , (123) 
and the matrix Q, given by 



Q 



-d x + d x V -d y 



-d„ 



dyV 



dyV 
d X +d X V 



(124) 



dx sin 2 (6») dx 



Then the matrix Tt = Q T Q is diagonal, with the operator 
of Eq. (|121|) as diagonal elements. The construction of 
a chiral matrix H is therefore complete for this problem, 
too: H in fact has the form of a random Dirac Hamil- 
tonian, studied in Ref. |5fJ. While we have not been able 
to find a quantity for the two-dimensional system that 
plays the role of the partial energy in one dimension, we 
have formulated a generalization of the Burgers equation, 
which is sufficient to show that a solution to Eq. (|123fl ex- 
ists. The practical determination of the chiral potential, 
however, remains an open problem in this case. 



Introducing coordinates for deviations from the ground- 
state configuration, 6 = 9q +ipg and <j) = 0o+V>) we take 
as canonical variables tpg and ip^ sin (#0). The quadratic 
Hamiltonian can be written 7i — Q T Q, with Q given by 
Eq. (|lll|l in terms of the 2x2 matrix V, which can be 
computed from £ using 

\ cot(8 ) + ^4nr^d 2 £ ^Lnded^S - 2^Hd j,£ 

This formalism would provide the starting point for a 
treatment of the one-dimensional random field Heisen- 
berg model similar to that presented for the XY model 
in Sec. 



VIII. RELATION TO EXPERIMENT 

The range of physical systems to which the ideas we 
have set out may apply is potentially very wide. The 
essential requirements are: quenched disorder; a ground 
state that can be considered classically; and excitations 
that can be treated as weakly interacting normal modes. 
. These requirements may be met for vibrational modes, 
either in glasses, or in randomly pinned phases with bro- 
ken translational symmetry such as charge density wave 
states. They may also be satisfied for magnetic excita- 
tions in disordered ferromagnets or antiferromagnets, or 
in spin glasses, either with or without significant mag- 
netic anisotropy. The most specific feature for which 
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one would like experimental evidence is probably the fre- 
quency dependence p(ui) oc oj 4 of the density for excita- 
tions that are not Goldstone modes. Alternatively, for 
excitations that are Goldstone modes, interest focusses 
on the frequency dependence of the mean free path at 
small oj, or, in the case of Heisenberg antiferromagnets 
and spin glasses which are quasi one-dimensional, the 
low-frequency form of the density of states. 

Before summarising the experimental situation, it is 
useful to outline ways in which excitations in a disor- 
dered system may fall outside the framework we have 
used. Large quantum fluctuations are the obvious route 
to different physics, and may be important in at least two 
ways. First, it can happen that a ground state is very 
far from being a classical configuration dressed with small 
zero-point fluctuations: random singlet phases^ in disor- 
dered quantum spin chains constitute a well-studied ex- 
ample. Second, it may be that even low-lying quantized 
excitations are unlike weakly interacting bosonic modes. 
To have a concrete example, consider the single mode 
problem for an anharmonic oscillator with potential U (q) , 
as reviewed in Sec. lIVBl Typically, this potential energy 
will have local minima separated by barriers from the ab- 
solute minimum, and the importance of quantum motion 
through or over the barrier will depend on the size of the 
mass m appearing in Eq. (|38l) . Thus, for a given U(q), in 
the semi-classical limit the low-lying quantum states will 
be close to harmonic oscillator levels with the classical 
normal mode frequency, while if quantum fluctuations 
are large, tunneling through barriers may hybridise lev- 
els associated with different classical minima, generating 
two-level systems, or the low-lying levels may be deter- 
mined by the form of U (q) in regions too far from its ab- 
solute minimum for the classical normal mode frequency 
to be relevant. In the ensemble, this crossover takes place 
as a function of frequency, as discussed in Ref. 37. As a 
result, one expects harmonic excitations with a density 
p(u) oc oj 4 at higher frequencies, and two-level systems 
with a constant density at lower frequencies. 

Turning to experiments, somewhat surprisingly, the 
best evidence that we are aware of for an excitation den- 
sity with an oj dependence is from studies of vibrational 
excitations in glasses. Here, inelastic neutron scattering 
and other measurements^ indicate an excess density of 
harmonic modes, compared to what is expected from De- 
bye theory and the measured speed of sound. Analysis of 
the temperature dependence of the heat capacity^ sepa- 
rates two contributions that are additional to the Debye 
value. One is approximately linear in temperature T, 
dominates at low T, and is attributed to two-level sys- 
tems. The other varies as T 5 , dominates at higher T, 
and is attributed to harmonic modes with the stated fre- 
quency dependence for their density. From a theoretical 
viewpoint, such excitations differ in an important way 
from the one we have discussed, because they coexist 
with propagating phonon modes. The consequences of 
coupling between the two sets of modes remain to be 
investigated, although coupling between localised anhar- 



monic vibrational modes and phonons has been studied 
in Ref. 

A context in which localised vibrational modes are ex- 
pected without any coexisting propagating phonons is 
provided by pinned charge density waves, represented in 
one dimension by the model studied in Sec. lVll Indeed, 
it was in this framework that Aleiner and Ruzin2£ and 
Fogler— argued for a density of harmonic proportional to 
oj 4 , with a crossover at low frequency^ to a constant den- 
sity of two-level systems. As reviewed by Fogler— , exist- 
ing measurements of frequency-dependent conductivity^ 
do not show the response expected from such harmonic 
modes, possibly because the low-temperature limit is not 
accessed. 

Studies of spin waves in disordered magnetic systems 
present opportunities to examine both Goldstone and 
non-Goldstone modes. In particular, inelastic neutron 
scattering measurements of spin dynamics in a dilute 
near-Heisenberg antifcrromagnet- ■ show the expected 
broadening in wavevector of excitations with increas- 
ing frequency. In contrast, magnetic neutron scattering 
measurements on amorphous magnetic alloys in which 
there is local magnetic anisotropy££ find modes which 
are broad in wavevector at all frequencies. In this case 
the density of excitations is approximately constant in 
frequency over the measured range. From our results, we 
expect a decrease in this density at low frequency, and 
a more detailed examination of low-frequency behavior 
would be of considerable interest. 



IX. CONCLUDING REMARKS 

Since we have investigated a variety of different direc- 
tions in this paper, it is perhaps useful to close with a 
short summary of our main points. 

Considering the general quadratic Hamiltonian for 
bosonic excitations, Eq. i[7)l. we have noted a formal sim- 
ilarity between it and fcrmionic random matrix Hamilto- 
nians from the additional symmetry classes (e.g. Eq. JSJ)). 
We have argued that this has limited direct consequences, 
because stability of the ground state requires correla- 
tions between matrix elements of the bosonic Hamilto- 
nian which invalidate a random matrix approach. In- 
stead we have shown that there is a useful mapping to 
an auxiliary problem with the structure of the chiral sym- 
metry class, which we have set out explicitly for a range 
of models. 

Examples of bosonic excitations in disordered media 
separate into those that are Goldstone modes, and those 
that are not. For the former, we have reviewed estab- 
lished results, which demonstrate that low frequency ex- 
citations decouple from disorder except for some systems 
in and below a marginal dimension d c = 2. For exci- 
tations that are not Goldstone modes, we have under- 
lined the way in which disorder itself generates low fre- 
quency excitations, and the universal form expected for 
their density, p{oj) oc oj 4 . Taking as an example the one- 
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dimensional XY model in a random field, we have used 
our techniques in detailed analytical and numerical cal- 
culations, obtaining results that illustrate this behavior 
of p(uj). 

A further application of the techniques we have devel- 
oped here is to Mattis models. These model spin glasses 
lack frustration, and their statistical mechanical proper- 
ties are equivalent under a gauge transformation to those 
of ferromagnets. They have a ground state spin configu- 
rations which are known for each disorder configuration, 
but excitations are nevertheless affected in a non-trivial 
way by disorder^ Because the ground state is known 
explicitly, the magnons in Mattis glasses are much eas- 
ier to study than in a real spin glass. Being Goldstone 
modes, these magnons fall into the same category as 



those in weakly disordered antiferromagnets, discussed 
in this paper. In particular, their critical dimension is 
d c = 2. In future work«2i by one of the present authors 
and A. Altland, the localization and transport properties 
of magnons in two and three dimensional Mattis glasses 
will be explored. 
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